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prehensive regional oceanographic data synthesis and integration effort for the U.S. Pacific Northwest 
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world-wide operational data centers were evaluated and summarized to provide scientific information and 
insights on system variability and processes affecting the PNW and, in particular, the OCNMS. The synoptic- 
scale viewpoint provided by satellite data made it possible to combine multiple observations on physical, 
biological, and sediment characteristics to fill existing data gaps, and document oceanic events for the entire 
region. The result is a unique representation of PNW ocean climate and behavior, an oceanographic base- 
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Program, and other Federal and State agencies, universities, and tribal staff. 

The characterization consists of three complementary components: a text report; a suite of quantita- 
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EXECUTIVE SUMMARY 

This report presents the results of a two-year investigation and summary of oceanographic satellite data ob- 
tained from multiple operational data providers and sources, spanning years of operational data collection. 
Long-term summaries of Sea Surface Temperature (SST) and SST fronts, Sea Surface Height Anomalies 
(SSHA), surface currents, ocean color chlorophyll and turbidity, and winds are provided. 

Overall Conclusions 

Merged satellite oceanographic data revealed information on: (1) seasonal cycles and timing of transition pe- 
riods; (2) linkages between seasonal effects (warming and cooling), upwelling processes and transport; and 
(3) nutrient/sediment sources, sinks, and physical limiting factors controlling surface response for Olympic 
Coast marine environments. These data and information can be used for building relevant hind cast models, 
ecological forecasts, and regional environmental indices (e.g. upwelling, climate, "hot spot") on biological 
distribution and/or response in the PNW. 

Key Findings 

(1) Patterns of SST and SSHA reveal linkages between seasonal heating effects and dominant physical pro- 
cesses affecting the sanctuary including upwelling, downwelling, and transport. 

SST patterns reveal seasonal maximums in August and minimums in January, with spring and fall transi- 
tions in May and November. Juan de Fuca Strait outflow results in suppressed SST during summer. Up- 
welling influences on SST are more dramatic inside the sanctuary versus points south of the sanctuary. 
SSHA patterns reflect vertical changes from mean sea level and reveal seasonal maximums in January 
and minimums in July, coincident with the peak of the downwelling and upwelling seasons, respectively. 
Maximum transition periods occur in March and November. Heights respond to coastal wind forcing, up- 
welling-induced (salinity and temperature driven) effects and freshwater influence inside the sanctuary. 
• Altimetry-derived radar estimates of SSHA could be used in combination with thermal infrared SST to 
isolate upwelling activity and biological response for PNW marine environments. 

(2) Patterns of ocean color chlorophyll and turbidity reveal insights on nutrient/sediment sources, timing of 
bloom response, and wind-driven processes affecting the sanctuary 

Chlorophyll patterns show spatially distinct pulses during the spring bloom period (April-June) just north 
of and inside the sanctuary, followed by strong upwelling-induced pulses inside the sanctuary from July 
through mid-October. 

Turbidity maximums occur during the fall-winter period inside the sanctuary; minimums occur in associa- 
tion with the Juan de Fuca Strait outflow in summer. 

Columbia River plume influence on ocean color is most evident in southern portions of the study area, 
with seasonal precipitation significantly correlated with ocean color turbidity inside the sanctuary. 

(3) Co-variability of surface water properties reveals subseasonal processes and modes of influence that 
limit chlorophyll production during the downwelling/upwelling cycles. 

Correlation patterns among precipitation, discharge, and ocean color in January suggest a strong coastal 
influence on surface sediment response characteristics inside the sanctuary. 

Correlation patterns among SST, SSHA, and ocean color in May suggest that increased regional heating, 
physical circulation and mixing processes limit chlorophyll production most dramatically. 
Correlation patterns among SST, SSHA, and ocean color in July suggest that increased light and tem- 
perature, along with upwelling, drive production inside the sanctuary. 

Correlation patterns in September suggest that coastal processes (Juan de Fuca outflow and Ekman 
transport) lead to fall chlorophyll maximums inside the sanctuary. 
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(4) Climatological anomalies of surface water properties reveal varied levels of response to environmental 
perturbations. 

Broad-scale climate patterns such as El Nino Southern Oscillation (ENSO) promote varying levels of 
surface anomalies inside the sanctuary, including extreme positive and negative SST, SSHA, chlorophyll, 
and turbidity. 

Regional and local OCNMS SST and SSHA anomalies trend positive with ENSO warm phases and nega- 
tive with cool phases. 

Regional and local OCNMS chlorophyll and turbidity anomalies vary in both strength and direction with 
documented ENSO and non-ENSO (neutral) periods. 
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1. INTRODUCTION AND BACKGROUND 

In January 2007, the Center for Coastal Monitoring and Assessment (CCMA), Office of National Ma- 
rine Sanctuaries (ONMS), and Federal/Regional Integrated Ocean Observation System (IOOS) Associations 
initiated discussions on the importance of broad-scale use of satellite-based oceanographic data to support 
the Olympic Coast National Marine Sanctuary (OCNMS). Satellite data requirements were identified to: (1) 
fill existing gaps in baseline environmental data, (2) enhance the understanding of ecological patterns and 
processes for the Olympic Coast region, and (3) provide information and insights to support the Olympic 
Coast Sanctuary Management Plan Review process (http://olympiccoast.noaa.gov/). To support these re- 
quirements, CCMA conducted a two-year investigation and summary of ocean climate, character, and color 
for the Pacific Northwest (PNW) region through satellite remote sensing. 

This oceanographic summary represents a unique regional synthesis, interpretation, and character- 
ization of multiple long-term satellite observational datasets made available through the National Aeronautics 
and Space Administration (NASA), National Oceanic and Atmospheric Administration (NOAA), and other 
worldwide operational data centers (see data products and sources). Observational monitoring emphasizes 
real-time data collection and archiving, product generation, and distribution of information vital for improved 
understanding of physical and biological characteristics in a regional context. A key benefit of sustained data 
collection is the potential to provide a clear synoptic-scale vision of important temporal and spatial oceanic 
events, useful for developing environmental indicators of perturbation and for identifying long-term trends 
within U.S. (and its associated territories') marine protected areas, including sanctuaries. 

The objectives of the characterization were to: 

(1) Identify, obtain, and perform regional syntheses on pertinent long-term satellite data to characterize the 
physical and oceanographic conditions within the PNW, and specifically the OCNMS; and 

(2) Summarize surface patterns, features, and underlying processes governing oceanographic variability and 
climate as depicted through satellite remote sensing for consideration in National Marine Sanctuary biodi- 
versity/conservation evaluations, management plan reviews, ecosystem assessments, and other decision- 
support initiatives. 

Herein, seasonal and interannual variability (and co-variability) of surface water properties and fea- 
tures were characterized with reference to dominant physical processes for the region. Seasonal surface 
patterns and changes over a typical annual cycle were analyzed and documented. Specific spatial-temporal 
linkages among water properties were analyzed to provide insights on ocean behavior and timing of sea- 
sonal transitions. Ocean temperature fronts and chlorophyll variability measures were used to reveal spatial- 
temporal aspects of phytoplankton (biological) activity for the region. In addition, oceanographic anomalies 
of key surface properties were documented as potential indicators of environmental perturbation, and for 
evaluating future trends for the sanctuary. 

This study complements prior oceanographic investigations through the use of temporally and spa- 
tially robust remote sensing datasets. A characterization of oceanographic baseline conditions for the OC- 
NMS is among the most important components in management decision support because it sets the context 
for describing the distribution of biological resources and for defining management alternatives. Ultimately, 
the goal is to link long-term synoptic imagery with direct observational and interpolated data to improve the 
understanding of environmental variability and perturbations affecting the sanctuary, and to monitor changes 
occurring in the sanctuary. 

The OCNMS was established in 1 994 in recognition of its nationally significant and unique marine en- 
vironment. The sanctuary resides in a dynamic oceanographic setting with broad scale (> 500 km) variability 
dominated by seasonal upwelling and downwelling. As a result, the PNW and OCNMS support rich marine 
mammal and seabird faunas, abundant populations of kelp and intertidal algae, and productive fish and 
invertebrate communities (ONMS 2008, Sydeman and Ellitot 2008). Due to the dynamic nature and harsh 
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oceanographic conditions along the PNW, collecting physical and biological monitoring data poses a signifi- 
cant challenge when using shipboard in situ and/or moored buoy monitoring. The satellite-derived proxies 
developed as fully georeferenced map layers were used to fill existing data gaps in baseline environmental 
information, and to provide key insights on local variability for the PNW region. 

1.1 Regional Setting and Study Area 

The study area lies in the northern portion of the Oregonian biogeographic province that extends from 
Pt. Conception, California, to Cape Flattery, Washington (Airame et al. 2003). The province is characterized 
by a narrow continental shelf, mountainous shoreline, steep rocky headlands, sandy pocket beaches with 
sea stack islands, many small and a few large rivers, and small estuaries with barrier islands. The province 
is also noted as exhibiting the greatest volume of upwelling in North America. The study region resides within 
the California Current System (CCS) and represents one of North America's most productive marine ecosys- 
tems. 




Vancouver Island 
j£5a CANADA 




Bamfield, BC 



The OCNMS, located in the northern portion of the study area off the coast of Washington, encom- 
passes approximately 8,572 km 2 (3,310 mi 2 ) of coastal and oceanic waters, or approximately 10% of the 
study area for this assessment (Figure 1). The sanctuary boundary extends from Koitlah Pt. due north to the 
U.S. -Canada border seaward to the 100 fm isobath. Following the 100 fm (-183 m) isobath southward, the 
boundary extends approximately 145 km, cutting across the heads of Nitnat, Juan de Fuca, and Quinault 
Canyons to a point due west of the Copalis River. The shoreward boundary is the mean lower low water line 
when adjacent to Indian reservations and state and county lands. When adjacent to federally managed lands, 
the boundary extends to the mean higher 126 o W 125 o W 124 o W 

high water line. 

1.2 Regional Oceanographic Summary 



The dominant physical broad-scale 
processes influencing PNW oceanography 
have been summarized by Hickey (1979, 
1 989, 1 998), Hickey and Banas (2003), Strub 
et al. (1987), and Strub and James (1988, 
2002). The PNW region is characterized by 
distinct patterns in oceanographic circula- 
tion, winter storms, flow-topography interac- 
tions, and land-sea influences. Large- scale 
processes are the predominant influence 
on seasonal upwelling-downwelling fluctua- 
tions, resulting in a highly dynamic oceano- 
graphic environment. 

PNW oceanography is dominated by 
the CCS, a broad, meandering system of cur- 
rents with strong interannual, seasonal and 
weekly variability (Hickey 1998). The CCS 
forms the southerly flowing segment of the 
greater North Pacific Gyre and extends from 
British Columbia, Canada to Baja California, 
Mexico (Figure 2). The dominant scales of 
circulation over much of the CCS are set up 
by strong alongshore winds, bottom topog- 
raphy, and shelf influences (Halliwell and Al- 
len 1987). Because of these characteristics, 
coastal trapped waves are efficiently gener- 




Figurel. Physiographic features off Washington coast and within 
Olympic Coast NMS. 
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ated and propagate long distances along the continen- 
tal margins of much of western North America (Hickey 
and Banas 2003). 

The CCS includes the slow, southward-moving 
California Current, the wintertime, northward-flowing 
Davidson Current, California Undercurrent, and other 
shelf and slope currents with primarily shorter-than- 
seasonal time scales. The southward-flowing California 
Coastal Jet Current and northward Davidson current 
are commonly observed on the continental shelf during 
the summer and winter, respectively (Figure 

3). The California Current flows southward year round 
offshore from the shelf break to a distance of 1 ,000 km 
from the coast (Hickey 1979, 1998). 



Figure 2. General schematic of the California Current System 
and North Pacific Gyre. Adapted from NOAA/NCCOS. 



Topographic features of the continental shelf in- 
teract with and influence flows and circulation in the 
region. The surface outflow patterns from the Juan de 
Fuca Strait, and interactions with bottom topography, 
affect chemical and biological processes, including 
transport mechanisms for nutrients, larvae, plankton, 

and even toxicants (domoic acid) associated with the Juan de Fuca eddy (Trainer et al. 2002, Trainer and 

Suddleson 2005, Trainer et al. 2009). 

The PNW includes one major river plume (the Columbia River), and several smaller estuaries and 
plumes, particularly evident south of the Juan de Fuca Strait. The Columbia River plume is highly dynamic 
with variation in direction, depth, and width in response to wind strength and direction (Thomas and Weath- 
erbee, 2006). The plume flows northward over the shelf and slope in the fall and winter. In the spring and 
summer it flows southward, well offshore of the Oregon shelf. 



1.3 Seasonal Cycles 

Seasonal oceanographic cycles influencing the PNW have been described in detail by Hickey and 
Banas (2003) and are defined by the dominant circulation patterns in the PNW region. Generally, oceano- 
graphic seasons can be characterized as (1) Winter (November-mid-February), (2) Spring (late February- 
March), (3) Early summer (April-June), (4) Late summer-early fall (July-September), and (5) Fall (October- 
November), including spring and fall periods of transition. 

Coastal upwelling is the dominant process controlling Olympic Coast water property variability over 
several-day to seasonal time scales. Upwelling along the Olympic Coast typically exemplifies colder, denser, 
saltier, and nutrient-rich waters that are being forced upward adjacent to the coast, typified by southward 
winds driving an offshore surface water movement known as Ekman transport. This occurs along the entire 
U.S. West Coast, but the strength and duration of upwelling increases to the south. With the exception of re- 
gions affected by the Columbia plume, stratification in the CCS is remarkably similar at most locations and is 
largely controlled by seasonal changes in large-scale advection (transport) and upwelling of water masses. 
Increases in chlorophyll production and biomass are evident during the heightened upwelling season. The 
region experiences periods of wind relaxation when seas are intermittently calm and winds are variable. 

Coastal downwelling generally occurs during the fall and winter, where a marked reversal from the 
summer pattern is evident along the Olympic Coast. The strong Davidson Current affecting the region is 
driven not only by the large-scale north-south pressure gradient, but also by strong winds from the south as- 
sociated with storms entering the region. Wind stress combining with the Coriolis force is especially strong 
in winter, driving an onshore Ekman transport of water (Thomas, Strub, and Brickley 2003). Downwelled wa- 
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ter along the coast is typically warmer, less saline, and 
nutrient reduced. The region can undergo diminishing 
cross-shelf and alongshore currents, causing a decrease 
in onshore Ekman transport at times during the winter sea- 
son (Allen and Newberger, 1996; Austin and Barth, 2002). 

For this report, seasonal oceanographic cycles and 
transition periods will be characterized by surface patterns, 
and changes in response to the dominant circulation driv- 
ers affecting the Olympic Coast region. Seasonal cycles 
will be documented to better understand expected condi- 
tions of each variable and to provide insights on dominant 
processes influencing variable response. General lags as- 
sociated with surface characteristics in response to ocean 
processes also will be documented. 
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1.4 Interannual Cycles 

Interannual and decadal-scale variability of the 
CCS involves complex atmospheric interconnections as- 
sociated with El Nino Southern Oscillation (ENSO), the 
Pacific Decadal Oscillation (PDO) and the North Pacific 
Gyre Oscillation (NPGO) (Halpin et al. 2004, Di Lorenzo et 
al. 2008). These cycles promote significant long-term vari- 
ability in the region that can yield extreme oceanographic 
anomalies on interannual to decadal time scales (Schwing 
et al. 1996). ENSO conditions are a result of interannual 
changes in sea level pressure between the eastern and 
western hemispheres of the tropical Pacific (Conlan and 
Service 2000). PDO is a long-term climatic pattern that is 
expressed by altering sea surface temperatures, surface 
winds, and sea level pressure (Mantua and Hare 2002). 
PDO shifts generally occur every 20-30 years and have 
occurred five times in the last century. (Airame et al. 2003). 
The PDO can be described as an extensive El Nino-like 
pattern of Pacific climate variability with warm (positive) 
and cool (negative) phases. The most recent shifts from 
warm/cool and cool/warm occurred in 1998 and 2003, 
(Airame et al. 2003; Peterson and Schwing 2003). During 
warm PDO, cold surface water temperature anomalies are 
observed in the western and central North Pacific Ocean, 
while the eastern Pacific exhibits above-average tempera- 
tures. The opposite patterns occur during cool phases. The NPGO is a climate pattern driven by broad-scale 
variations in ocean upwelling and transport processes. The NPGO shows strong interannual patterns of 
association with key physical-biological ocean variables, including salinity, nutrients, and chlorophyll a (Di 
Lorenzo etal. 2008). 
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Figure 3. Generalized surface flow of the California Cur- 
rent System in the PNW during typical (a) summer and (b) 
winter conditions. Source: NOAA. 



For this report, interannual patterns of surface variable response will be presented as climatological 
anomalies (deviations) to highlight spatial-temporal characteristics of those deviations from expected condi- 
tions. Variable response will also be visually compared with documented shifts in ENSO. 
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2. DATA AND METHODS 

2.1 Data Synthesis 

Remote sensing data provide key information to understand ocean environments by characterizing 
the distribution and variability of oceanic surface features. These include ocean fronts, chlorophyll production 
and transition zones, cold and warm core eddies, convergence zones, and river plumes. Optical observa- 
tions of physical, biological, and sediment characteristics provide valuable quantitative information on local- 
ized transport processes, meso-scale circulation features, and even large-scale water mass characteristics. 
Radar altimetry estimates of sea surface heights and roughness readily depict ocean surface characteristics 
including mean sea level, marine currents, winds, and waves. The remote sensing data presented herein 
were synthesized based on 8+ years of operational coverage, and include weekly to monthly averaged, and 
climatological anomaly data, for satellite-derived sea surface temperature (SST), sea surface height anomaly 
(SSHA), geostrophic currents, ocean color products, and winds. In addition to remote sensing data, gridded 
monthly land precipitation and Columbia River discharge data were used for comparison with remote sensing 
data, where appropriate. A complete product index for the data used throughout this report can be found in 
Tablel. 

Sea Surface Temperature (SST) 

The Coral Reef Temperature Anomaly Database (CoRTAD) SST estimates were obtained from the 
NOAA/NASAAdvanced Very High Resolution Radiometer (AVHRR) Oceans Pathfinder Version 5.0 Program, 
developed by the University of Miami's Rosenstiel School of Marine and Atmospheric Science (RSMAS) and 
the NOAA National Oceanographic Data Center (NODC). The Pathfinder V 5.0 Program measures global 
cloud cover, sea surface temperature, and ice, snow, and vegetation cover and characteristics. Weekly 
(7-day) and monthly averaged data for both the ascending pass (daytime) and descending pass (nighttime) 
are available on equal-angle grids of 8,192 elements in longitude and 4,096 in latitude (4 km resolution). 
Gap-filled weekly averaged 4km data were used to achieve the greatest possible spatial coverage over a 
7-day cycle. The CoRTAD SST "baseline" temporal resolution encompasses 1 ,096 weeks for the period Jan- 
uary 1985 to December 2005. CoRTAD data were then subset to include the entire West Coast, and weekly 
data were then binned to monthly averages. CoRTAD SST data series are useful for a range of applications 
including climatological (1985-2005) baseline and trend determination, ocean frontal boundary delineation, 
thermal stress detection in coral reefs, and climate indicator development (Halpern et al. 2008, Bruno et al. 
2007). 

SST Fronts 

SST fronts were detected using CoRTAD data as described above, but also from NOAA/NESDIS/ 
STAR Office holdings of the NOAA Geostationary Operational Environmental Satellites (GOES-10/11) SST 
and SST fronts data. SST fronts typically occur at the boundaries between various surface water types char- 
acterized by different physical, biological and optical properties. Frontal probabilities were calculated from 
daily-averaged SST fields based on 24-hourly observations. SST fields were processed at NOAA/STAR 
and CCMA by applying an edge-detection algorithm to identify SST fronts. The probability of detecting fronts 
(PDF) was calculated based on number of times a particular pixel qualifies as a front divided by the number 
of clear pixels (following Castelao et al. 2006). 

Sea Surface Height Anomaly (SSHA) 

SSHA data from altimetry reflect vertical changes from mean sea level and were obtained from the 
Archiving, Validation and Interpretation of Satellite Data in Oceanography (AVISO) Program. AVISO delayed 
time products of sea surface height generated from merged Topex/Poseidon (T/P), Jason-1/2, ERS-1/2 and 
ENVISAT missions. Delayed time (DT) Reference series products were used to provide accurate, homoge- 
neous measurements across time. Weekly (7-day) and monthly averaged data are available on equal-angle 
grids of 1080 elements in longitude and 720 in latitude (referred to as 1/4° grid size). In the derivation of 
SSHA, an estimate of mean sea level (referred to as the mean profile from T/P data for the period January 
1993 to December 1999) was applied to develop one merged vertical anomaly map combining measure- 
ments from available altimeter missions. SSHA data used in this report encompassed 729 weeks from Oc- 
tober 1992 to September 2006. Monthly estimates of vertical sea surface height from mean sea level were 
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Table 1. Remote sensing data products and source specifications for the Olympic Coast National Marine Sanctuary 
Oceanographic Characterization. 



Product Type 


Data Source 


Frame 


Resolution 


Summary 


Units 


Pathfinder SST 
and SST Anomaly 
(CoRTAD) 


NOAA/NESDIS/NODC 


1985-2005 


4km 


Weekly/ 
Monthly 


°C 


SST gradients 
(CoRTAD) 


NOAA/NESDIS/NODC 


1985-2005 


4km 


Weekly/ 
Monthly 


(°C/12km) 


GOES-1 0/11 SST 
and SST fronts 


NESDIS/NODC/STAR 


2000-2007 


4km 


Daily/ 
Monthly 


°C 


SeaWiFS Ocean 
Color Chlorophyll 
and Anomalies 


NASA 


1997-2007 


1km 


Daily/ 
Monthly 


MgL-1, 
Steradian -1 


SeaWiFS Ocean 
Color Turbidity; 
Remote Sensing 
Reflectance at 
670nm) 


NASA 


1997-2007 


1km 


Daily/ 
Monthly 


Steradian -1 


Jason-1, 
Topex/Poseidon, 
ERS-1/2ENVISAT 
Sea Surface 
Height Anomalies 


AVISO, 

SSALTO/DUACS & 
CNES 


1992-2006 


1/4° grids 


Weekly/ 
Monthly 


cm 


Jason-1, 
Topex/Poseidon, 
ERS-1/2ENVISAT 
Geostrophic 
Surface Currents 


AVISO, 

SSALTO/DUACS & 
CNES 


1992-2006 


1/3° grids 


Weekly/ 
Monthly 


cm/s, 

degrees from 
north 


QuikSCAT Sea 
Surface Winds 


Remote Sensing 
Systems 


1999-2007 


25km 


Monthly 


m/s, degrees 
from north 


Land Precipitation 
Totals and 
Anomalies 


GPCC 


1986-2006 


1° grids 


Monthly 


mm 


Columbia River 
Streamflow Totals 


USGS 


1986-2006 


watershed 


Monthly 


cfs 



obtained following a simple bin averaging technique. Merged SSHAdata series provide information on sea 
level trends, density-driven changes, geostrophic current intensity and direction, and subsurface layers, i.e., 
mixed and thermocline depths (Strub and James 2002, Wilson and Coles 2004, Gordon and Giulivi 2004). 

Geostrophic Currents 

Estimates of geostrophic currents were obtained from merged Jason-1, Topex/Poseidon, ERS-1/2, 
and ENVISAT missions provided by the AVISO Program. Delayed time Mean Sea Level Anomaly (MSLA) 
products of geostrophic currents were used. Using oceanographic convention, weekly (7-day) Easterly-to- 
ward the East (U) and Northerly-toward the North (V) components, and monthly averaged U-V data, were 
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made available on a Mercator grid projection with a resolution of 1/3° grid size, encompassing 729 weeks 
from October 1992 to September 2006. Geostrophic currents result from the balance between Coriolis force 
and horizontal pressure gradients. Geostrophy implies that the slope of SSHA measured from satellites is di- 
rectly related to the pressure gradient at the sea surface, and thus, to a geostrophic velocity at the surface. 

Ocean Color 

The ocean color dataset from Sea-viewing Wide Field-of-View Sensor (SeaWiFS) provided a 10-year 
(September 1997 to October 2007), near daily set of chlorophyll and light attenuation data, i.e., turbidity, at 
1-km spatial resolution. SeaWiFS Ocean Color Level-2 archive products produced and distributed by the 
NASAGoddard Space Flight Center's Ocean Color Data Processing System (OCDPS) include the geophysi- 
cal values for chlorophyll a using the OC4v5 algorithm. Chlorophyll a is the dominant pigment in marine pho- 
tosynthetic organisms, and is referred to simply as chlorophyll within this report. Light availability, or turbidity, 
is estimated using the remote sensing reflectance at 670nm Rrs(670) in units of Steradian-1 , which is a proxy 
for light scattering due to inorganic sediments in the water. SeaWiFS data are useful for investigating pat- 
terns, trends, and episodic events of near-surface ocean phytoplankton variability over a multitude of spatial 
domains (Sackman et al. 2004, Thomas, Townsend, and Weatherbee 2003). 

Winds 

QuikSCATscatterometerwind products derived from Remote Sensing Systems, Inc. (RSS) provided 
a 7-year (July 1999 to September 2006), weekly set of wind magnitude and direction data at 25km spatial 
resolution. Winds were calculated using a microwave scatterometer to estimate sea surface roughness and 
radiometer measurements for rain flagging and sea ice detection. Weekly and monthly averaged data were 
derived from all data within the calendar month. Wind speeds are scalar averaged and directions are vector 
averaged. 

Land Precipitation and PNW Drainage Boundaries 

Regional precipitation data used for spatial comparisons were obtained via the Global Precipita- 
tion Climatology Centre (GPCC) real-time precipitation data product based on the complete GPCC monthly 
rainfall station database in 1° x 1° resolution (Adler et al. 2003). GPCC land precipitation monitoring data 
products were processed using area-averaged and time-integrated precipitation fields, based on surface rain 
gauge measurements collected nationally. The geographic framework (taken from NOAA/NOS/SPO 1999) 
depicts the boundaries of basins that would include the majority of runoff entering the shelf area of the Olym- 
pic Coast, including those from river systems and coastal drainages. Only portions of GPCC grid cells that fell 
within the coastal drainage basins were used. Monthly precipitation totals were combined as a mean value 
for the entire coastal basin. 

Climatological Summaries 

Seasonal to inter-annual shifts in oceanic and atmospheric conditions off the Olympic Coast can 
significantly influence the persistence, location, and spatial orientation of surface layer features and spatial 
patterns. The satellite image data series allowed for the extraction of weekly to monthly surface data and 
for the development of long-term (climatological) weekly to monthly data summaries of average (expected) 
conditions for each variable presented. In addition, monthly temporal deviation (climatological anomalies) 
and variability (standard deviation or index) terms were calculated for each variable (V). The general formula 
used in the determination of expected conditions and anomalies follows: 



Vi\rn,y) = V\(rn,y) - V\{m) 

where, V l represents the monthly average (or median) for month m of year^, V\ is the long-term (grand) 
mean for month m of all years, and V* is the climatological anomaly for month m of year y. This calculation 
was performed on all overlapping pixels in the monthly and weekly image databases to create image cli- 
matologies for each variable. These climatological images and anomalies will be used to summarize water 
properties for the upwelling, downwelling, and transition periods, and over interannual cycles. 

Assessing variability of surface water properties from satellite provides added spatial and temporal 
context into the physical processes that limit or enhance phytoplankton production at the surface and at depth. 
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Figure 4. Transects used for seasonal satellite oceanographic 
data extractions. 



To address chlorophyll variability, the range (spread) 
around the monthly and climatological monthly me- 
dian (50 th percentile) values were calculated using the 
difference between the 75 th and 25 th quantiles for each 
overlapping pixel in each month of the calendar year. 
Quantiling is used to eliminate questionable values in 
the SeaWiFS dataset. Resulting maps show the spa- 
tial distribution of months with the highest variability, 
revealing chlorophyll "patches" with shape, size, and 
position occurring at different months of the year. 

2.2 Visual Representation of Data 

The development of consistent daily and 
weekly time series of multiple image datasets allowed 
for visual representation and interpretation across 
variables. Since monthly averaged images are es- 
sentially snapshots of expected or average conditions 
over time, and are not representative of within month 
variability, extractions of weekly data were performed 
to help explain temporal patterns of variability impor- 
tant to ocean climate and character within and adja- 
cent to the sanctuary. Data were extracted along three transect lines in the study area, following a general 
longitudinal (E-W), cross-shelf, and longshore (N-S) orientation (Figure 4). The longitudinal E-W transect was 
used to compare patterns of variability inside the sanctuary versus the open ocean environment; the diagonal 
cross-shelf transect was used to compare shelf patterns influenced by Juan de Fuca Strait outflow and eddy 
circulation versus shelf slope and open ocean; and the latitudinal N-S transect was used to compare long- 
shore patterns inside the sanctuary versus areas influenced by estuarine and river plumes to the south. 

Long-term weekly averaged SST (Figure 5), SSHA, chlorophyll, and turbidity data were extracted 
along the three transect lines to evaluate spatial-temporal patterns of variability and linkages with comple- 
mentary satellite datasets. Vector quantities of SSHA-derived geostrophic currents extracted into northerly 
(toward the north) and easterly (toward the east) components were temporally averaged and plotted as 
conventional stick plots, depicting magnitude and direction of surface currents (Figure 5, bottom). Weekly (in- 
stantaneous) and anomaly data were extracted using time-series color plots (Hovmoller diagrams) and line 
plots to reveal interannual characteristics of variables (Figure 6). A Hovmoller diagram maps a scalar quantity 
(e.g., SST in °C) over a span of latitude or longitude through time. The scalar quantities were extracted in the 
longitudinal, cross-shelf, and longshore directions to depict spatial-temporal variability and episodic (anoma- 
lous) events for the period. 

2.3 Quantitative Representation of Data 

Statistical measures of association between remote sensing data types were used to highlight physi- 
cal and biological coupling, and to infer the nature of oceanographic processes along the Olympic coast. Be- 
cause these large data sources represent co-located and broad-ranging surfaces (e.g., maps of sea surface 
temperature, chlorophyll concentration), correlation was used to assess the magnitude and direction (posi- 
tive/negative) of association among oceanographic variables. Because many of these variables were not 
normally distributed, Spearman's nonparametric distribution-free rank correlation was used (Lehmann and 
D'Abrera 1998). This technique provides a measure of monotone (linear) association, and is often used when 
the distribution of the data make Pearson's correlation coefficient undesirable or misleading. Correlation 
coefficients (p - rho) were computed and mapped among monthly mean/median data, and for climatological 
anomaly data. Correlation matrices describing spatial and temporal associations of all variables, as well as 
inferred mechanisms influencing variable associations were assessed and documented. 
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Data Disclaimer 

All remote sensing data, particularly near the coast, should be treated with caution. Near-shore al- 
timetry data that are extrapolated from offshore data typically result in poor sea surface height estimates and 
subsequently poor geostrophic velocity calculations (Saraceno et al. 2008). Ocean color chlorophyll signals 
may be biased in some cases due to the presence of color-dissolved organic materials (CDOM) and particu- 
late matter in coastal regions (Gregg and Casey 2004). Turbidity reflectance fields also can be affected by 
the presence of algal pigments in the water column due to absorption of chlorophyll in the red wavelengths 
(DairOlmo and Gitelson 2005). 
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Figure 5. Example plot showing long-term weekly averaged (expected) SST conditions along a North-South transect line for 
Olympic Coast shelf waters (red line) and Washington Coast shelf waters (blue line), and SSHA-derived surface currents depicting 
magnitude and direction along the southern portion of the transect. Time period for CoRTAD SST: 1985-2005; AVISO geostrophic 
currents: 1992-2006. 
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Figure 6. Example time-series plot (Hovmoller) of weekly gap-filled CoRTAD SST (C°) revealing SST along a North-South transect 
line for Olympic Coast shelf waters (red line) and Washington Coast shelf waters (blue line). Time period: January1995-December 
2005. 
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3. RESULTS 

In this section, satellite ocean climate data were summarized to reveal seasonal and interannual vari- 
ability and transition inside and outside of sanctuary boundaries. Transition is defined as a typical seasonal 
shift in surface patterns for each variable presented. Highlighted will be the unique spatial and temporal char- 
acteristics and differences across variables for the three transects analyzed. Next, the correlation patterns 
across select remote sensing data types will be evaluated to help identify the scales of influences (local to 
regional, ephemeral to semi-permanent) on OCNMS water property response. Visual tools including Hov- 
moller diagrams, correlation (rho) maps, and statistical plots will be presented to illustrate the significance of 
these relationships. Finally, interrelationships across all data types will be utilized to reveal important limiting 
factors and modes of influence on surface layer response inside and outside of sanctuary boundaries. 

3.1 Seasonal Patterns of SST, SSHA, and Currents 

SST data were examined to reveal seasonal characteristics for the PNW region, addressing spatial 
and temporal variability inside and outside the sanctuary. Seasonal SST patterns (Figure 7) reveal cool win- 
tertime temperatures (purple- blue), indicative of the downwelling season, followed by gradual spring transi- 
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Figure 7. CoRTAD SST monthly climatologies (1985-2005) showing typical SST distribution during winter (December-February), 
spring (March-May), summer (June-September), and fall (October-November) periods. 
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tion in March and April. May reveals increased surface heating influence from the south, especially along 
the coast (cyan-green). Surface heating in summer results in dramatically increased SST mean fields from 
early June through September, temporally associated with the upwelling season. SST averages range from 
8° C to near 16° C in a 4-month interval. A fall transition in SST typically occurs in October and November. 
Climatological SST maximums occur during August and climatological minimums occur in February (Figure 
7). 

Seasonal SSHA data were summarized to evaluate coastal perturbations indicative of wind-forc- 
ing, density changes from heating/cooling, and salinity (upwelling effects). Because areas surrounding the 
OCNMS are significantly affected by strong seasonal variation in coastal currents, dynamic month-to-month 
transitions in SSHA are expected (Figure 8). Seasonal SSHA patterns show positive (orange-red) nearshore 
heights during the downwelling season, 8-1 Ocm above the vertical mean. Circulation changes during late 
winter to early spring result in significant SSHA transition from February to April (Figure 8). Spring months are 
characterized by mainly neutral (cyan) patterns evident in April. During summer upwelling, a sign reversal (+ 
to -) from winter occurs inside the sanctuary, where heights are negative (blue), 10-1 2cm below the vertical 
mean. A gradual fall transition occurs in October and November. SSHA maximums inside the sanctuary oc- 
cur during January and minimums occur primarily in June/July, with a secondary peak in September. 
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Figure 8. AVISO SSHA monthly climatologies (1 992-2006) showing typical SSHA distribution during the winter (December-Feb- 
ruary), spring (March- May), summer (June-September), and fall (October-November) periods. 
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East-west summary 

Hovmoller diagrams of expected 
weekly SST and SSHA reveal the dramat- 
ic spring SST increases (cyan to green) 
up until a mid-late summer maximum 
(orange-red), followed by the dramatic 
cooling effect in November (Figure 9). A 
significant transition occurs starting in May 
(cyan to green approximates a transition of 
1 .5 C°; Figure 9). A similar SST transition 
maximum occurs in November. Warmer 
waters penetrating northward in spring re- 
veal earlier transition in near-shore areas 
than in the open ocean (Figure 9, top). Fall 
transition is more gradual nearshore than 
in offshore areas (Figure 9, bottom). 

SSHA signal strength and sign 
(+/-) are highly dynamic near shore versus 
offshore. A wintertime maximum (yellow- 
orange; Figure 9) near shore occurs dur- 
ing January with a peak of approximately 
7cm above mean sea level. In offshore 
areas, SSHA seasonal variability is more 
gradual, trending towards a September- 
October maximum. SSHA minimum inside 
the sanctuary typically occurs during mid- 
late September, coincident with increased 
upwelling at that time. Spring surface 
height transition occurs before SST tran- 
sition with maximum sea-level transition 
in March (Figure 9). Fall sea-level tran- 
sition occurs in November. SSHA-de- 
rived surface geostrophic currents reveal 
north-northwestward flow in winter and 
southward flow in summer/fall. Strongest 
currents occur in late September and Jan- 
uary. 

Cross-shelf summary 

Hovmoller diagrams illustrate gen- 
erally lower SST associated with the Juan 
de Fuca Strait outflow and the eddy fea- 
ture during summer (Figure 10). Mean 
SST along the transect reveal summer 
temperatures to be 14.6°C (red line; fig- 
ure 10), approximately 2°C colder than in 
areas of the continental slope (blue line). 
Lower SST in shelf waters forms an SST 
frontal boundary in early fall, generally as- 
sociated with the 200 m isobath. 

Strong near-shore variations in 
SSHA signals are apparent along the 
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Figure 9. Climatological summary of weekly averaged (expected) conditions 
of SST and SSHA during a typical annual cycle along an East-West longitudinal 
(47°N) transect line for the continental shelf (<200m isobath; red line) and con- 
tinental slope, open ocean (>200m isobath; blue line). SSHA-derived general 
surface currents are depicted from extracted data along the shelf break of the 
study area (box). 
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cross-shelf transect, where summer min- 
imums can exceed 10cm below mean 
sea level in areas outside of the Juan de 
Fuca Strait (Figure 10). Of note are the 
pronounced maximum heights in Janu- 
ary, especially in nearshore areas, where 
heights can exceed 10cm above the 
vertical mean. Decreased salinities and 
coastal impingement are speculated as 
the predominant drivers to elevated sea 
heights. Spring sea-level transition typi- 
cally occurs two months earlier than SST 
transition (Figure 10; boxes). Fall sea- 
level transition is typically one month af- 
ter SST transtion. Maximum differences 
between near-shore and offshore heights 
are slightly exaggerated in January and 
September (Figure 10), coincident with 
stronger surface currents during these 
months. Geostrophic currents along 
the continental slope (Figure 10) reveal 
north-northwestward flow in winter, and 
southward flow in summer/fall. 

Alongshore summary 

Hovmoller diagrams reveal no- 
table differences in SST and SSHA inside 
the sanctuary, with cooler waters persis- 
tent in June-July just south of Bamfield, 
BC, near the mouth of the Juan de Fuca 
Strait (Figure 11). Fall SST minimums 
occur in September just south of Cape 
Flattery. Wintertime SST differences be- 
tween northern and southern sections of 
the transect are negligible, with a Febru- 
ary minimum of 8.5°C. 

Longshore variations in SSHA 
patterns are evident inside the sanctuary, 
with a maximum in January and mini- 
mums throughout the upwelling season. 
A mainly neutral SSHA pattern is evident 
near the mouth of the Columbia River. Of 
note are the dramatic differences between 
SST minimum and SSHA maximum in 
January. Fluctuations in sea level are 
influenced by changes in temperature, 
salinity, and winds. The heightened sea- 
sonal SSHA patterns in winter suggest a 
strong salinity and wind-driven influence, 
especially since the colder wintertime 
seawater is less dense than it is during 
summer. 
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Figure 10. Climatological summary of weekly averaged (expected) conditions of 
SST and SSHA during a typical annual cycle along a cross-shelf transect line 
for the continental shelf (<200m isobath; red line) and continental slope, open 
ocean (>200m isobath; blue line). SSHA-derived general surface currents are 
depicted from extracted data along the shelf break of the study area (box). 
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Figure 1 1 . Climatological summary of weekly averaged (expected) conditions of SST and SSHA during a typical annual cycle along 
a North-South transect line for Olympic Coast shelf waters (red line) and Washington Coast shelf waters (blue line). SSHA-derived 
general surface currents are depicted from extracted data along the shelf break of the study area (box). 

Summary of seasonal processes influencing SST and SSHA 

Patterns of SST and SSHA for the Olympic Coast reflect responses caused by a combination of phys- 
ical factors. During a strong upwelling period in September 2001 , regional SST minimums corresponded with 
sea surface height minimums (Figure 1 2). The suppressed SSHA summer signals are indicative of upwelling 
pulses injecting cooler dense waters to the surface, creating lower sea surface heights (Figure 12; top row). 
Lower SST is generally shoreward of maximum equatorward flow (Figure 1 2). During the downwelling period 
of January 2004 (Figure 12; bottom row), northward wind-induced coastal impingement, Juan de Fuca out- 
flow, and coastal/estuarine influence were likely to have resulted in positive SSHA patterns inside the sanctu- 
ary. A more thorough investigation of processes influencing surface temperatures and sea level fluctuations 
will be presented in section 3.4 and 3.5. 

Depictions of both SST and SSHA, derived from infrared and radar altimetry estimates, respectively, 
reveal dominant seasonally induced processes affecting the region, i.e., wind-driven upwelling and down- 
welling. Since radar altimetry is not cloud-dependent, subtleties in SSHA variability are captured in the AVISO 
data, which are sometimes lost in the SST signals due to cloud obscurement. SSHA appears to be less re- 
sponsive to summer air/sea heating and more responsive to upwelling-induced mixing and wind influences. 
In summary, these variables could be used in combination to isolate upwelling response for the region. 
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Figure 12. Monthly mean SST distribution, SSHA, and SSHA-derived surface current 
velocities for September 2001 and January 2004, during the height of the upwelling and 
downwelling periods, respectively. 
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3.2 Seasonal Patterns of Ocean Color Chlorophyll, Turbidity and Winds 

Ocean color data were used to reveal seasonal characteristics for the PNW region, addressing spa- 
tial and temporal variability inside and outside of sanctuary boundaries. The climatological monthly mean 
images shown illustrate the dynamic nature of the spatial distribution of monthly and seasonal chlorophyll 
(Figure 13). Lower chlorophyll (cooler colors) is indicative of the downwelling season (Figure 13; top row) , 
most apparent during December through February, followed by a gradual spring transition starting in March. 
In May, coastal chlorophyll concentrations north of the sanctuary increase dramatically, exhibiting a peak in 
the 6.0- 1 0.0 |jg/L range during the start of the upwelling season. During August to October, a primary peak 
occurs inside the sanctuary, when highest chlorophyll typically appears for the season. A fall transition typi- 
cally occurs in November, when chlorophyll production significantly diminishes, falling to lows of around 1.5 
|jg/Lfrom November through January (Figure 13; bottom row). 
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Figure 13. SeaWiFS monthly climatologies (1997-2007) showing typical chlorophyll distribution during the winter (December-Febru- 
ary), spring (March), early summer (April-September), late summer-early fall (July-September) and fall (October- November) peri- 
ods. 
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Ocean color turbidity, as measured by remote sensing reflectance at 670nm (Rrs670) is a proxy for 
light availability and particulate (inorganic) loading in the upper water column (Figure 14). This surface effect 
can be caused by runoff, wind-generated resuspension of sediments, or upwelling currents. Turbidity values 
are generally very low for the PNW; high values are relatively uncommon and spatially limited. Regions of 
persistently high values are associated with river plumes and the sediment generated by coastal runoff west 
of the Olympic Mountains. The highest median turbidity values are the lighter colors (greens to yellows), and 
lower values are depicted by cooler colors decreasing from greens to blues (Figure 14). The patterns for 
light availability are divided into two main seasons — wet and dry. Coastal areas south of the Juan de Fuca 
Strait experience their highs in the winter wet season when rainfall is at its highest, while areas to the west of 
the Strait, along the continental shelf, experience a primary peak in summer (Figure 14), during the dry sea- 
son. Values for near-shore areas in the sanctuary are highest from December through March, with a spring 
transition evident in April, falling to a seasonal low in May. The fall months of September and October are 
characterized by transition, marking the end of the dry season and start of the wet season. 

East-West summary 

Hovmoller plots of SeaWiFS data extractions along the E-W transect reveal expected chlorophyll 
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Figure 14. SeaWiFS monthly climatologies (1 997-2007) showing typical Rrs670 value distribution during the winter wet season (De- 
cember-February), spring (March-April), summer dry season (May-September), and fall (October- November) periods. 
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concentrations and Rrs670 values for 
the climatological period of January- 
December (Figure 15). Values repre- 
sent expected conditions based on the 
10-year climatology. A climatological 
minimum near shore occurs from No- 
vember through early February, where 
the width of coastal chlorophyll is re- 
strained to extreme coastal areas and 
concentration peaks at ~5 |jg/L. Slight 
increases in concentration are evident 
in late February, with significantly stron- 
ger peaks beginning in late July through 
early October inside the sanctuary (Fig- 
ure 15, red line). In July, the width of 
coastal chlorophyll widens significantly. 
Climatological maximum near shore oc- 
curs during late September-early Octo- 
ber, when maximum concentrations can 
reach 8 |jg/L. Of note are the chloro- 
phyll responses offshore in mid-October 
extending into the open ocean environ- 
ment, where chlorophyll can increase 
from 0.1 |jg/L to 0.5 |jg/L during a two- 
week period (Figure 15, blue line). The 
onshore-offshore chlorophyll enhance- 
ment ends fairly abruptly in November. 

The differences in remotely 
sensed signals of turbidity versus chlo- 
rophyll are evident in Figure 15. From 
November through March, a band of 
turbidity > 0.002 Steradian-1 can be 
seen during the wet season. Between 
the end of April to mid October, turbid- 
ity near the coast reduces by an order 
of magnitude to -0.0005 Steradian -1, 
with a minimum ocurring in May. Mini- 
mal turbidity increases are evident near 
the shelf break in summer, with coastal 
increases building in late October and 
peaking in December. 

Winds show strong seasonal 
patterns with east-south-eastward di- 
rections in spring and early summer, 
south-eastward directions in mid-sum- 
mer to early fall, and northward winds in 
winter (Figure 15). 
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Figure 15. Climatological summary of weekly averaged (expected) conditions of 
SeaWiFS chlorophyll, turbidity and QuikSCAT winds during a typical annual cy- 
cle along an East-West longitudinal (47°N) transect line for the continental shelf 
(<200m isobath, red line) and continental slope, open ocean (>200m isobath, blue 
line). Winds were extracted along the transect line near the shelf break (box) and 
vectors plotted to reveal direction in oceanographic convention. 
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Cross-shelf summary 

Through pixel extractions along 
the cross-shelf transect, Hovmoller 
diagrams illustrate spring chlorophyll 
blooms as well as the influences of the 
Juan de Fuca outflow in spring and sum- 
mer. Chlorophyll along the coast is near 
1 |jg/L from December through March. 
Starting in early April, chlorophyll in- 
creases of up to 6 |jg/L are observed 
(Figure 16). The chlorophyll maximums 
indicative of the spring bloom period can 
be seen (Figure 16, red line). Nutrient 
buildup in the cross-shelf area is likely 
influenced more strongly by spring river 
runoff out of the Juan de Fuca strait at 
this time. Of note are the fall peaks evi- 
dent off the shelf break in October prior 
to the fall transition period (Figure 16). 

A spatially homogeneous area 
of more turbid water persists just inside 
of the shelf break from May to Septem- 
ber. Slight increases in Rrs670 returns 
during upwelling could be caused by 
resuspended sediments churned up by 
the strong upwelling currents. Of note 
is the near-shore area around the Juan 
de Fuca Strait outflow, which exhibits 
decreased turbidity starting in June and 
ending in October (Figure 16). 

Winds show stronger magni- 
tudes offshore, with eastward directions 
shifting southeastward in summer and 
early fall (Figure 16). During fall transi- 
tion, a fairly abrupt shift from eastward 
to northward occurs in November and 
December. 

Alongshore summary 

Longshore variability in chlo- 
rophyll signals along the N-S transect 
line highlight the distinct spring bloom 
periods to the north of the sanctuary, 
the summer/fall peaks inside the sanc- 
tuary, and the Columbia River outflow 
influence to the south (Figure 17). Dur- 
ing fall transition (November) and win- 
ter downwelling, chlorophyll decreases 
significantly near shore. A distinct tran- 
sition occurs in late February along the 
mid-transect, where general increases 
are seen over a two-week period, and 
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Figure 16. Climatological summary of weekly averaged (expected) conditions of 
SeaWiFS chlorophyll, turbidity and QuikSCAT winds during a typical annual cycle 
along a cross-shelf transect line for the continental shelf (<200m isobath, red line) 
and continental slope, open ocean (>200m isobath, blue line). Winds were extract- 
ed along the transect line off the shelf break (box) and vectors plotted to reveal 
direction in oceanographic convention. 
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then a flattening occurs until late March (Figure 17). By early April, distinct blooms consistently >6 |jg/L are 
evident in northern sections (Figure 17, red line), lasting until July. By August, the northern enhancement 
edge dissipates, with more southerly enhancement persisting through mid October, with maximums ranging 
from 8-10 |jg/L inside sanctuary boundaries. Of note is the implied influence of nutrients from the Columbia 
River plume affecting the nearshore, remotely sensed chlorophyll signals in the immediate vicinity of the 
plume (Figure 1 7). Also of note is the apparent break in the elongated signal band and decreased chlorophyll 
of <0.5 |jg/L predominantly inside the sanctuary in November and December. 

Turbidity values in the alongshore direction reveal influences of the wet season in fall and winter 
(Figure 17). From the second half of December to the end of March, the magnitude on the southern end of 
the transect is -0.002 1/sr. While this magnitude will remain constant until August, the northward extent of 
coastal influence from the Columbia River varies. The northward extent increases from December through 
February, and then recedes in March through May. The northbound extent dissipates until the last week of 
August, and then drops an order of magnitude to 0.0005 1/sr. The strong precipitation influence west of the 
Olympic Mountains may be attributed to the increased Rrs670 signal in November and December, where 
land-sea runoff may play a significant role in elevated Rrs670 distribution. 

Winds reveal three prevailing patterns, with slightly lighter winds in spring and early summer, stonger 
upwelling favorable winds southeastward in mid-summer, and northward winds in fall and winter (Figure 17). 
Fall wind shifts to the north typically occur during October after the fall primary chlorophyll peak. 
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Figure 17. Climatological summary of weekly averaged (expected) conditions of SeaWiFS chlorophyll, turbidity and QuikSCAT 
winds during a typical annual cycle along a North-South transect line for Olympic Coast shelf waters (red line) and Washington 
Coast shelf waters (blue line). Winds were extracted along the transect line along the shelf (box) and vectors plotted to reveal direc- 
tion in oceanographic convention. 
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The Juan de Fuca Strait outflow and circulation 
eddy has been shown to be a significant source of nu- 
trients for the region (Hickey and Banas 2003, Trainer 
et al. 2002). Using SeaWiFS Rrs670 values represent- 
ing typical low conditions (less turbid) during summer, a 
narrow band of less turbid waters is evident in areas just 
outside and to the south of Cape Flattery (Figure 18, 
purple). These lower Rrs670 values could represent a 
proxy to nutrient source for this area. 

Summary of ocean color response characteristics 

Ocean color chlorophyll and turbidity estimates 
provide information on algal production, light availabil- 
ity, and river/estuarine plume influence inside and out- 
side of the sanctuary. Chlorophyll conditions during a 
strong upwelling period in May 2002 revealed spatially 
distinct pulses outside of the Juan de Fuca Strait with 
slight increases in Rrs670 values (Figure 19). Strong fall 
upwelling-induced pulses were evident during Septem- 
ber 2006 inside sanctuary boundaries, with fairly stable 
Rrs670 values (Figure 19). Chlorophyll maximum typi- 
cally occurs along the Olympic Coast in late Septem- 
ber and early October. Although Rrs670 values are low 
year-round, values generally exhibit a cyclical temporal 
variability, with late fall-winter peaks generally occurring 
during the precipitation maximum for the region. Colum- 
bia River plume influence on ocean color is most evident 
in southern sections of the study area. Juan 
de Fuca outflow in summer induces clearer 
water conditions inside the sanctuary. May-2002 
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Figure 18. SeaWiFS-derived climatological estimates of 
Rrs670 20th quantile values during July revealing decreased 
turbidity associated with Juan de Fuca outflow inside sanctu- 
ary boundaries. Note bifurcation near the mouth of the Strait. 
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Figure 19. Monthly mean SeaWiFS chlorophyll and turbidity (Rrs670 val- 
ues) for May 2002, and September 2006 during the spring bloom period 
and enhanced fall upwelling periods, respectively. 



An Oceanographic Characterization of the Olympic Coast National Marine Sanctuary 



3.3 OCNMS Seasonal Cycles and Transition Periods 

Through analysis of surface patterns of chlorophyll, currents, and winds, five seasonal cycles were 
documented for the Olympic Coast region. Seasonal cycles and transition periods were noted for sanctuary 
waters only, and may or may not be indicative for other areas outside the sanctuary. 

(1) Winter (November through mid-February) 

(2) Spring transition (late February through March) 

(3) Spring/early summer bloom period (April through June) 

(4) Late summer/early fall transport and upwelling period (July through early October) 

(5) Fall transition (October through November) 

(1) Winter (November through mid-February) — This period is dominated by strong winds from the south, on- 
shore Ekman transport, northward mean flow, and downwelling. The present study depicts surface patterns 
consistent with previous investigation summarized by Hickey and Banas (2003) and Hickey (pers. comm). 
Chlorophyll patterns are lowest for the year in December (present study). Currents and winds reveal pre- 
dominantly north-northwest patterns during the height of the downwelling cycle (Figure 20). The region can 
also undergo diminishing cross-shelf and alongshore currents, causing a reduction in onshore transport at 
times during the winter season (Allen and Newberger 1996; Austin and Barth 2002). 

(2) Spring transition (late February through March) — This period reveals variable winds and currents, with 
generally mean northward flow. Surface pattern climatologies reveal general shifts in currents, winds and 
chlorophyll production during this time (Figure 20). Upwelling of deep nutrient-rich waters begins after the 
transition and continues through the fall (Landry et al. 1989). An "early spring" period is evident in a February 
chlorophyll pulse depicted inside the sanctuary, and could be representative of a transition of fair, upwelling- 
favorable wind conditions, followed by increased volatility and storminess. A somewhat distinct 2-3 week 
westward wind shift is evident (Figure 20). 

(3) Spring/early summer bloom period (April through June) — During this period, upwelling continues, bringing 
deep waters to the surface near the coast in some areas. Surface flows over much of the shelf and slope are 
generally southeastward, especially later in the period (Figure 20). Spring chlorophyll enhancement occurs 
as upwelled nutrients reach the surface near the coast and in the Strait of Juan de Fuca and its eddy (Mac- 
Fadyen et al., 2005; Hickey and Banas, 2003). The period is associated with increasing SST regionally. 

(4) Late summer/early fall transport and upwelling period (July through early October) — This period is domi- 
nated by southeastward mean surface flow (Hickey 1989, present study), Ekman-induced cross-shelf trans- 
ports, and upwelling. Coastal upwelling and transport are the dominant processes controlling Olympic Coast 
water property variability during this time (Lentz 1992, Kirincich et al. 2005). Increases in chlorophyll pro- 
duction are evident during this period. Both currents and winds are generally lighter during this time (Figure 
20). 

(5) Fall Transition (October through November) — During the fall transition, winds shift to more downwelling- 
favorable (northward), and mean monthly surface currents begin to transition more northward with corre- 
sponding changes in water-column stratification (Hickey 1989). The upwelled waters relax downward offshore 
across the shelf and nutrient supply thus diminishes. A fairly abrupt decrease in chlorophyll concentration 
occurs during this time. 
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Figure 20. Diagrams of seasonal breakouts and periods of transition for the Olympic Coast Sanctuary. Average (expected) 
weekly SeaWiFS median chlorophyll data extracted along a North-South transect line, AVISO geostrophic currents (box), and 
QuikSCAT winds (circle) depict weekly patterns over the calendar year. Breakouts are characterized by surface chlorophyll 
patterns and circulation drivers for the Olympic Coast region. 
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3.4 Seasonal Patterns of Oceanic Fronts, Gradients, and Ocean Color Variability 

Ocean fronts are boundary features typically associated with intensified motion, sharp gradients, 
decreased stability, and increased convergence at the ocean surface. Efforts to characterize fronts using 
satellite-derived edge detection methods are well documented (Castelao et al. 2006, Breaker et al. 2005). In 
PNW marine environments, fronts can be placed into three general classes: (1) small-scale fronts associated 
with estuarine/river discharge, (2) upwelling fronts associated with strong thermal and salinity gradients, and 
(3) open ocean fronts caused by wind-forcing. Fronts in the PNW are most apparent during and slightly after 
the heightened upwelling season. 

In OCNMS ocean environments, fronts initially develop in May (Figure 21), coincident with SST un- 
dergoing the greatest temporal transition. Coastal upwelling intensity increases dramatically in June, and 
fronts strengthen during this time. By July, frontal intensity is centered inside the shelf break near the Juan 
de Fuca eddy circulation, and by September, persistent fronts set up near the shelf break. The most intense 
offshore fronts occur in October and November. Frontal activity inside the sanctuary decreases significantly 
during the November transition, reaching a minimum between February and April (Figure 21). Small-scale 
fronts are also associated with estuarine/river discharge. 
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Figure 21. GOES climatological (2000-2007) SST monthly frontal probabilities showing long-term frontal persistence during 
the winter (December-February), spring (March-May), summer (June-September), and fall (October-November) periods. 



An Oceanographic Characterization of the Olympic Coast National Marine Sanctuary 



Ocean color chlorophyll variability estimates the level of change of chlorophyll concentration over a 
monthly interval. Assessing variability from satellite data provides added spatial context into ocean phyto- 
plankton activity and physical processes that limit or enhance biological production at the surface. Using dif- 
ferences between the 75 th and 25 th quantiles, resulting maps show the spatial distribution of months with the 
highest variability, revealing distinct and sometimes pronounced areas of enhancement. 

At the height of the "turbid season" in January, chlorophyll variability is lowest for the year. A distin- 
guishable area of increased chlorophyll variability is evident in the February image (Figure 22). April patterns 
of heightened chlorophyll variability exhibit a narrow band of heightened variability inside the sanctuary (Fig- 
ure 22). Variability then decreases in the sanctuary in late summer, when upwelling pulses are speculated 
as the dominant process of influence. In September, during the peak of seasonal chlorophyll production, vari- 
ability also increased and was more widespread, mainly confined outside of sanctuary boundaries (Figure 
22). 
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Figure 22. SeaWiFS- derived climatological (1 997-2007) estimates of chlorophyll variability during the winter (December-Feb- 
ruary), spring (March), early summer (April-September), late summer-early fall (July-September) and fall (October-November) 
periods. Variability was determined using the difference of the 75th and 25th quantiles for the distribution of chlorophyll and 
turbidity values, respectively. 
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Ocean color variability 
speaks to enhanced transition in 
biotic response. It is assumed that 
chlorophyll response can be mea- 
sured not only through algal abun- 
dance (using chlorophyll a as a 
proxy), but through temporal change 
in abundance, which may be more 
useful for coupling with other pa- 
rameters. Figure 23 depicts these 
heightened chlorophyll transition 
periods effectively. Specifically high- 
lighted are the changes associated 
with the time periods of late March 
(Figure 23a), depicting a calm pe- 
riod, and mid-April (Figure 23b), 
showing significant enhancement in 
chlorophyll variability. 
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Figure 23. Comparison of climatological chlorophyll variability estimates during a) 
late March, versus b) mid-April. Note the significant enhancement in variability 
inside the sanctuary. 



Juan de Fuca Eddy 

The Juan de Fuca eddy is a semi-permanent cyclonic feature located off the coast of Vancouver 
Island and the Juan de Fuca Strait. The extreme diversity in habitats and zones of convergence of ocean 
currents accounts for the rich biodiversity in the eddy. 

The counterclockwise cold-core eddy is situated west of the Strait of Juan de Fuca, southwest of 
Bamfield, BC. The eddy is a dominant feature of circulation around the OCNMS with a diameter ranging from 
50-1 00km, and is visible in summertime satellite imagery in seasonal minimum low SST and, generally, mini- 
mums and maximums in chlorophyll. The seasonal eddy is a result of the interaction between effluent from 
the Strait, southward wind-driven upwelling processes, and underlying topography — a spur of the Juan de 
Fuca submarine canyon (Hickey and Banas 2003). The eddy forms in spring and declines in fall (Freeland 
and Denman 1982). The eddy was defined by using density profiles (MacFadyen et al. 2005). Decreases 
in SST in the center of circulation are most apparent in July, when climatological SST gradient magnitudes 
are most distinct. In days prior to and during Eddy SST enhancement, light to moderate NW winds prevail in 
near-shore zones, indicative of fair weather conditions. 

A well-documented edge detection algorithm (Canny 1986) was used to delineate SST frontal gradi- 
ent edges associated with eddy boundaries using both the CoRTAD and GOES long-term data series. In ad- 
dition to gradients, frontal orientation was also computed to provide a reasonable proxy of circulation around 
the center of circulation. Climatological gradients and along-front directions reveal implied circulation around 
the eddy during climatological June and July (Figure 24), indicating a slight westward shift in eddy orientation 
from June to July in the long-term data series. Spatial patterns of frontal directions, climatological median 
chlorophyll, and eddy location are shown in Figure 25, with implied circulation around the center depicted 
through along-front direction. 

Summary of SST fronts and chlorophyll variability 

SST fronts for the region show distinct patterns indicative of the Juan de Fuca eddy circulation, up- 
welling fronts, offshore wind forced events, and the Columbia River plume effect. Ocean color chlorophyll 
variability estimates the level of change of chlorophyll over a temporal interval. Heightened monthly chloro- 
phyll variability is typically oriented along the eastern edges of SST frontal boundaries, hinting at increased 
convergence and enhanced biological activity. Westward shifts in eddy circulation are evident in the long- 
term SST data series. Spatial patterns of near-shore and offshore frontal boundaries and chlorophyll infer 
zones of increased upwelling, and possible habitat niches for biological assemblages. 
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Figure 24. CoRTAD SST climatological (1985-2005) monthly gradients and along-front di- 
rections for June and July. 
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Figure 25. CoRTAD SST climatological (1985-2005) monthly along-front directions and 
climatological monthly median chlorophyll for June and July. 
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3.5 OCEAN CLIMATE DATA SUMMARY 

In this section, previously shown Hovmoller diagrams of remotely sensed SST, SSHA, chlorophyll, 
turbidity, currents, and winds are presented in combination to visually depict average (expected) conditions 
for the climatological period covering January to December. Extracted data along the three transect lines 
were compared across the four variables, and inferences were made on the dominant factors influencing 
surface response. 

East-West Summary 

Suppressed near-shore SST values are evident starting in early June and ending in October, indica- 
tive of seasonal upwelling processes (Figure 26). Suppressed SSHA values inside the sanctuary (Figure 
26, cyan-blue) are likely the result of density-driven (temperature/salinity) processes related to seasonal up- 
welling. In winter, heightened SSHA values along the coast are likely caused by precipitation/runoff, the Co- 
lumbia River, and coastal impingement (wind forcing). Spring chlorophyll blooms typically start in mid-April. 
Chlorophyll pulses are much more pronounced in late summer and early fall, typically ending in mid-October 
(Figure 26). Wind-induced effects can lead to SST fronts in offshore areas with the most persistent frontal 
features setting up in October and November, coincident with increased chlorophyll production offshore. 
Slight influence of localized precipitation and sediment resuspension can be seen in expected turbidity condi- 
tions (Figure 26). Of note are the strong northwestward-dominated mean currents and winds in January and 
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Figure 26. Diagrams of averaged (expected) daily-weekly AVISO SSHA, CoRTAD SST, SeaWiFS chlorophyll and turbidity (Rrs670), 
AVISO geostrophic currents, and QuikSCAT winds along an East-West transect line over the calendar year. Plots reveal typical con- 
ditions of remotely sensed surface response characteristics for OCNMS waters. 
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southeastward currents in September (Figure 26). 

Cross-shelf Summary 

Cross-shelf patterns in ocean climate reveal distinct transitions between near-shore and offshore sur- 
face water property responses. Outflow from the Juan de Fuca Strait causes suppressed SST during summer 
(Figure 27). Negative SSHA in summer is likely caused by dips in temperature due to Juan de Fuca outflow in 
combination with upwelling processes (Figure 27). Spring chlorophyll blooms typically start in mid-April and 
diminish slightly by mid-summer. Wind-induced effects can lead to SST fronal development (frontogenesis) 
in offshore areas. Of note are the SST and SSHA spatial transition areas coincident with increased chloro- 
phyll production (and fronts) in October. Influence of the Juan de Fuca outflow can be seen in expected tur- 
bidity conditions promoting lower turbidities directly adjacent to the Strait (Figure 27). A summer/fall turbidity 
maximum is also noted near the shelf break (Figure 27). Of note are the southeastward-dominated currents 
during July, with a transition to more southward mean flow in August and September (Figure 27). 

Alongshore Summary 

Ocean climate patterns in the alongshore direction reveal seasonal warming, density- driven, and 
wind effects along the Olympic Coast. Outflow from the Juan de Fuca strait is evident, where SST differ- 
ences between extreme northern and southern sections typically reach 2-3°C during summer. Sea-level 
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Figure 27. Diagrams of averaged (expected) daily-weekly AVISO SSHA, CoRTAD SST, SeaWiFS chlorophyll and turbidity (Rrs670), 
AVISO geostrophic currents, and QuikSCAT winds along a cross-shelf transect line of each data series over the calendar year. Plots 
reveal typical conditions of remotely sensed surface response characteristics for OCNMS waters. 
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fluctuations between summer 
and winter are most exagger- 
ated along this transect. Wind- 
induced coastal impingement, 
Juan de Fuca outflow, and up- 
welling processes appear to 
have significant influence on 
SSHA signals in northern sec- 
tions, with decreased influence 
near the Columbia River mouth 
(Figure 28). From ocean color 
values, the evidence of spring 
bloom events is most apparent 
in northern sections of the study 
area from early April through 
mid-June (Figure 28). Near 
shore Ekman-induced pulses 
then dominate the region; this 
is most apparent during mid-Ju- 
ly through early October when 
chlrophyll reaches an annual 
maximum. Of note is the Co- 
lumbia River outflow influence 
on chlorophyll concentrations 
and turbidity in extreme south- 
ern sections, and heightened 
chlorophyll (low Rrs670) during 
a brief period in February (Fig- 
ure 28). February is typically 
a time of pre-transition, when 
frequent, severe cold-front 
passages could cause bottom 
resuspension and subsequent 
increased chlorophyll at the sur- 
face. Although Rrs670 signals 
are generally low to moderate 
year-round, a coastal runoff ef- 
fect is noted across northern 
and southern sections from 
November through December. 
General transition periods are 
depicted well with altimetry- 
derived geostrophic currents 
(Figure 28, bottom). Coastal 
currents become redirected 
from typical wintertime pat- 
terns (north-northwestward) to 
southeastward starting in early 
March through April. Summer- 
time southeastward flow re- 
mains fairly consistent through 
October until the November 
transition period. 
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Figure 28. Diagrams of averaged (expected) daily-weekly AVISO SSHA, CoRTAD SST, 
SeaWiFS chlorophyll and turbidity (Rrs670), AVISO geostrophic currents, and QuikSCAT 
winds along a North-South transect line of each data series over the calendar year. Plots 
reveal typical conditions of remotely sensed surface response characteristics for OCNMS 
waters. 
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3.6 Surface Property Relationships and Processes Affecting the OCNMS 

The relationship between surface layer properties provides information on bio-physical coupling, and 
local versus regional processes affecting OCNMS waters. For example, global studies investigating the co- 
variability of chlorophyll and SST have demonstrated the effects of large-scale climate cycles on heightened 
or suppressed signals over expansive ocean regions (Yoder and Kenelly 2003, McGillicuddy et al. 2001). 
In the Pacific Ocean basin, El Nino periods generally favor lower chlorophyll and higher SST, while La Nina 
years favor high chlorophyll and low SST. Localized spatial patterns of SST and chlorophyll have been linked 
to hemispheric climate cycles and have been used to identify feeding grounds for many commercially impor- 
tant species (Valavanis et al. 2004). In this section, monthly SST, SSHA, ocean color chlorophyll and turbid- 
ity, Columbia River discharge, and land-based precipitation data will be used to help define spatio-temporal 
correlation peaks during downwelling and upwelling periods. Multiple rho maps will be used to depict select 
spatial correlation patterns of interest in both local and regional ocean environments. 
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SST and SSHA Co-variability 

Through pair-wise correlation approaches using Spearman's rho (p), the association of monthly SST 
and SSHA image data values were tested. When all months were tested, SST/SSHA relationship was 
positive, indicating that thermal expansion properties have 
a strong influence on SST/SSHA relationship year-round, a. 
During upwelling in July (Figure 29a), slightly neutral to 
negative patterns are evident near the Juan de Fuca Strait 
outflow. The density-driven processes related to upwelling 
appear to reinforce a positive relationship in October (Fig- 
ure 29b). Seasonal heating (thermal expansion effects) is 
the likely contributor to positive SST/SSHA correlation re- 
gionally in both summer and winter. 

SST, SSHA, and Chlorophyll Co-variability 

The overall correlation pattern between SST and 
chlorophyll anomaly values reveal mostly negative p in 
open ocean waters with neutral patterns along the Olympic 
Coast. That is, cooler SST co-varies with increased chlo- 
rophyll regionally. Strength and direction of association 
between these two variables are generally moderate and 
negative when averaged over the entire upwelling period, 
but reveal unique spatial-temporal patterns within the up- 
welling season. 

At the start of upwelling season in May, distinct spa- 
tial patterns arise in SST/chlorophyll correlation fields (Fig- 
ure 30). May shows a distinct negative association along 
the Olympic Coast shelf, specific to areas along the Van- 
couver Island and Washington coastlines. This pattern rep- 
resents a time when solar radiation effects are in transition 
and upwelling pulses of nutrient-rich, cooler waters become 
more evident. Differences between offshore and nearshore 
pixels is evident when monthly mean SST and median chlo- 
rophyll values are extracted from selected clusters (boxes, 
Figure 30). Trend lines of fit are added to emphasize oppos- 
ing relationship in offshore (Figure 31a) and coastal (Figure 
31b) environments. The nearshore relationship (Figure 
31b) emphasizes low SST values (10.5°C) corresponding 
with a heightened chlorophyll signal up to 3|jg/L. SSHA 
and chlorophyll correlations (not shown) reveal similar pat- 
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Figure 29. Monthly SST vs. SSHA Spearman correlation 
coefficients (rho) showing strength and direction (+/-) of 
spatio-temporal relationship for the months of (a) July and 
(b) October. 
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terns in May although negative correlation fields are not 
as spatially homogeneous. 

July correlations along the Olympic Coast reveal 
typical summer conditions, with increases in light avail- 
ability and sea surface heating acting to limit chlorophyll 
production at the surface. Nutrient availability in July 
increases overall due to increased upwelling and other 
localized processes. An opposing pattern from May is 
shown in the July plots, with a negative trend offshore 
(Figure 32a) versus a positive trend near shore (Figure 
32b). Chlorophyll concentration reached 14|jg/L for ex- 
tracted pixels in the sample when temperatures were 
greater than 15°C. Rho was moderate to strongly posi- 
tive near-shore, and negative offshore (Figure 33). The 
overall pattern offshore is strongly negative. SSHA and 
chlorophyll rho maps (not shown) show a similar near- 
shore versus offshore relationship in strength and direc- 
tion. 

In September/October, SST/chlorophyll correla- 
tion patterns reveal negative association regionally (Fig- 
ure 34), with a small cluster of positive rho associated 
with Juan de Fuca outflow and eddy circulation inside the 
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Figure 30. Spearman correlation coefficients (rho) showing 
strength and direction (+/-) of spatio-temporal relationship 
between CoRTAD monthly SST and chlorophyll for May. 
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Figure 31. Scatterplot displaying the relationship between CoRTAD SST and SeaWIFS chlorophyll during climatological May for 
waters (a) outside and (b) inside the sanctuary. Values are extracted from boxes shown in Figure 30. 
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Figure 32. Scatterplot displaying the relationship between CoRTAD SST and SeaWIFS chlorophyll during climatological July for 
waters (a) outside and (b) inside the sanctuary. Values are extracted from boxes shown in Figure 33. 
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Figure 33. Spearman correlation coefficients (rho) showing 
strength and direction (+/-) of spatio-temporal relationship 
between CoRTAD monthly SST and chlorophyll for July. 
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sanctuary. This pattern would indicate that localized pro- 
cesses related to the Juan de Fuca eddy may have more 
influence on SST/chlorophyll relationship inside the sanc- 
tuary. Differentiating the dominant physical processes 
that limit chlorophyll and sediment characteristics will be 
addressed in the following sections. 

Nutrient/sediment Sources and Ocean Color C o - 

variability 

Outflow from the Columbia River plume, Juan de 
Fuca Strait, and other land-sea interactions have a signifi- 
cant effect on biological activity and sediment response 
in the OCNMS (Thomas and Weatherbee 2006). In this 
section, co-variability of chlorophyll, turbidity, precipita- 
tion, and discharge will be presented to highlight sources 
(and sinks) of nutrient and sediments impacting OCNMS 
marine environments. 

The Columbia River plume is highly dynamic with 
variation in direction, depth, and width of response (Hickey 
and Banas 2003). The Columbia River plume core flows 
northward on the inner shelf in fall and winter and south- 
ward in spring and summer, offshore of the Oregon shelf. 
Through spectral classification techniques, Thomas and 
Weatherbee (2006) depicted the Columbia River plume 
core to reach near sanctuary boundary waters in winter. 
Spearman's rank correlation of monthly Rrs(670) and Co- 
lumbia River discharge values exhibit a similar pattern 
during the winter period (Figure 35a). A distinct "plume" 
of high rho (Figure 35a, deep orange) is evident from the 
Columbia River mouth northwards along the coast. Of 
note is the negative association (deep blue) in coastal 
areas south of Cape Flattery (Figure 35a). 

Increased precipitation and runoff play significant 
roles in the transport and distribution of sediments and 
nutrients out of Olympic Coast river systems, estuarine 
plumes, and the Juan de Fuca Strait (Thomas and Weath- 
erbee 2006, Newton et al. 2003). A strong potential exists 
for precipitation-effect on sediment and nutrient charac- 
teristics along the OCNMS inner shelf. Through pair-wise 
approaches, a significant portion of turbidity in near-shore 
zones co-varies with coastal precipitation during the win- 
ter wet season (Figure 35b). These correlation patterns 
support the premise that seasonal runoff from areas to 
the north of the Columbia River plume contribute signifi- 
cantly to higher turbidity values in OCNMS waters, relative to the influence of the Columbia Plume itself. The 
results speak to other factors, such as more localized coastal precipitation, as substantial contributors to 
seasonal Rrs670 response. 

The association of precipitation and chlorophyll should differ based on the amount of light available 
for algal stimulation. It is speculated that negative precipitation/chlorophyll correlation in September (Figure 
36a) indicates that optimal light conditions (less precipitation) could yield more production at the surface. In 
October, the relationship is positive but limited to an area around the mouth of the Juan de Fuca Strait, includ- 
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Figure 34. Spearman correlation coefficients (rho) showing 
strength and direction (+/-) of spatio-temporal relationship 
between monthly SST and chlorophyll for September. 
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ing the sanctuary (Figure 36b). October is more indicative of increased precipitation, which could influence 
correlation patterns together with the possible influence of both the Strait and coastal runoff processes. 
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Figure 35. Spearman correlation coefficients (rho) showing strength and direction (+/-) of spatio- 
temporal relationship between SeaWiFS January turbidity and (a) USGS gauged Columbia River 
January discharge and (b) GPCC January precipitation. 



128°W 



126°W 



124°W 




Correlation Coefficient (rho) 



0.25 +0.25 +0.50 +l 



Figure 36. Spearman correlation coefficients (rho) showing strength and direction (+/-) of spatio- 
temporal relationship between SeaWiFS chlorophyll and GPCC precipitation for (a) September 
and (b) October. 
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Summary of intraseasonal (within season) processing influencing surface layer relationships 

To report on processes affecting surface layer responses inside and outside OCNMS waters, a cor- 
relation matrix revealing level and direction of association between remote sensing parameters and other 
physical-hydrometerological variables is presented (Table 2). Correlation patterns of surface water properties 
are spatially and temporally tied to upwelling/downwelling cycles (mixing processes), increased light avail- 
ability, and changes in land-based physical forcing mechanisms. In this section, the dominant mechanisms 
influencing surface layer characteristics inside the OCNMS will be interpreted from sign (+/-) and strength of 
selected relationships from Table 2, and from the climate summaries of expected conditions (Figures 26-28). 

For all months, a positive relationship is evident for SST and SSHA (row 1 , Table 2). Thermal expan- 
sion properties and upwelling reinforces positive SST/SSHA relationship. During the start of the upwelling 
season in May, SST/chlorophyll and SSHA/chlorophyll relationship is negative inside the sanctuary and posi- 
tive regionally (row 2, 3; Table 2). Nutrient-rich, colder water being forced to the surface during upwelling 
enhances negative SST/chlorophyll relationship. In July, the SST/chlorophyll relationship is highly influenced 
by localized heating effects, inducing a positive relationship. During September and October, SST/chloro- 
phyll relationship is negative at regional scales, suggesting that chlorophyll is heavily influenced by physical 
circulation (water mass characteristics) in addition to upwelling (mixing processes). The relationship between 
ocean color chlorophyll/turbidity and precipitation is positive during January inside the sanctuary (row 6, 7; 
Table 2). Columbia River discharge and turbidity reveals both a negative and positive relationship inside 
sanctuary boundaries (row 5, Table 2). 

Table 2. Levels of association between remote sensing parameters, Columbia River discharge, and land-based precipi- 
tation inside and outside OCNMS boundaries. Possible regional and local processes influencing spatial-temporal levels 
of association are added for reference. 
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Inference to dominant processes affecting relationships 

Thermal expansion and density driven processes reinforce positive relationship inside 
sanctuary and in open ocean environments. During upwelling, colder subsurface 
waters forced upwards induces suppressed SST and SSHA mainly confined to areas 
east of the shelf break. 



Physical circulation and mixing processes influence negative relationship during 
upwelling months. Chlorophyll is partially enhanced by seasonal heating effects in 
July inducing positive relationship. In May, a positive offshore relationship generally 
occurs due to warmer waters penetrating northward promoting algal stimulation. A 
split correlation pattern exists locally with Juan de Fuca outflow influencing negative 
relationship inside sanctuary boundaries. 

Same as above. Negative relationship during upwelling months possibly linked to 
geostrophic transport and advection processes. Additional data testing needed. 



Terrestrial sources (coastal and Columbia River) reinforce positive relationship in 
January. During dry season, relationship is mainly neutral to slightly positive. 
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Precipitation (coastal runoff) induces positive relationship in nearshore areas during 
January, October and November (not shown). 
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Precipitation (coastal runoff) induces positive relationship in nearshore areas during 
January and October. Juan de Fuca Strait outflow to the south in May and September 
possibly reinforces negative relationship during the dry season. 




Positive correlation (> 0.77 rho) 
Negative correlation (< - 0.77 rho) 
Not enough data or no relationship 
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Monthly differences in localized processes and limiting factors affecting surface layer relationships 
from remote sensing are visually depicted in map form (Figure 37). During January, surface layer response 
inside the OCNMS is heavily influenced by downwelling conditions (poleward currents, onshore Ekman 
transport), precipitation (coastal/estuarine runoff), and Columbia River discharge (Figure 37). The positive 
association of precipitation, discharge, chlorophyll and turbidity suggests that the turbid source is from estua- 
rine/coastal waters, mainly west of the Olympic Mountains. 

During May, surface layer responses inside the sanctuary are mainly influenced by strengthened 
upwelling and Juan de Fuca Strait outflow (Figure 37). The negative May correlation patterns of chloro- 
phyll and SST are indicative of these limiting factors, where colder, nutrient-rich waters are forced upwards 
through upwelling processes promoting algal stimulation at the surface. Negative association of precipitation 
and chlorophyll suggests that the nutrient source is not from coastal/estuarine runoff processes. Flow out of 
the Juan de Fuca Strait is less turbid during spring and summer, which would indirectly support the negative 
association between decreased precipitation (with increased light availability) and increased chlorophyll. Re- 
sults were inconclusive on the influence of the Columbia River during May, but southward circulation would 
support lesser affect inside the sanctuary. 

During July, surface layer responses are indicative of increased light and temperature, with available 
nutrients promoting chlorophyll production along the entire Olympic coast. Although upwelling persists in 
the region at promoting enhanced chlorophyll, SST appears to be more sensitive to seasonal heating than 
to mixing processes. A negative association between chlorophyll and precipitation suggests that blooms are 
induced by either upwelled water or Juan de Fuca Strait outflow (Figure 37). 
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Figure 37. Schematic representation of physical factors limiting nutrient availability and chlorophyll production 
during (a) winter (January), (b) spring (May), (c) mid summer (July), and (d) early fall (September) periods. 
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September chlorophyll peaks (annual maximums) inside the sanctuary are likely tied to localized pro- 
cesses, including Juan de Fuca outflow and eddy circulation (Figure 37). In fall, nutrient availability remains 
high, when Juan de Fuca outflow is circulated and recirculated, heated, and forced southward, causing 
increased production south of Cape Flattery (Figure 37). It is speculated that increases in Ekman transport 
leads to increased production at the surface in September/October. 

Through correlation results and analysis of expected surface characteristics for the OCNMS, the gov- 
erning physical limiting factors on surface layer responses were assessed and documented for the sanctuary 
(Table 3). Contribution levels for these limiting factors are based on spatial and temporal sign and strength of 
relation for the months of January, May, July, and September, and are relative to five broad-scale physical fac- 
tors: Columbia River, Coastal/Estuarine, Juan de Fuca Strait, Upwelling/Marine, and Light/Temperature. The 
Columbia River shows moderate influence in January, with lower influence for the remaining months (column 
1, Table 3). Coastal/Estuarine influence shows a high level of contribution on response characteristics for 
January and October, with moderate influence during July and September. Juan de Fuca Strait outflow and 
circulation shows significant influence during all upwelling months analyzed. Mixing processes associated 
with upwelling/marine influence is strong during spring, summer and fall. Light- and temperature-promoting 
conditions necessary for algal bloom activity are most pronounced in July, with moderate effect occurring 
in spring and fall. The information in Table 3 is intended to provide a process-driven conceptual framework 
on governing processes at low-frequency (monthly) time scales, but does not infer higher-frequency (daily- 
weekly) patterns of co-variability. 

In summary, using 12 calendar months and six different variables, nearly 100 correlation maps were 
produced and analyzed. Tables 2 and 3 condense these patterns and governing processes that are most ap- 
parent from visual interpretation and quantitative-based analysis. Many of the correlation patterns document- 
ed in this report support existing local knowledge of the physical and biological factors that limit production 
on seasonal time scales. However, the month-to-month differences in relation between variables reveal the 
importance of long-term satellite data series in further dissecting month-to-month variability in this dynamic 
ocean environment. The long-term relationship with remotely sensed bio-physical properties allows for more 
robust synoptic-scale interpretations and inferences to be made on estuarine-near shore interaction between 



Table 3. Relative level of contribution of physical limiting factors on surface layer response inside OCNMS boundaries. Deter- 
mination of contribution level based on analyis of expected conditions and intra-seasonal correlation patterns. 

Limiting Factors/Mode of Influence 





Columbia River Coastal/Estuarine Juan de Fuca Strait Upwelling/Marine Light and Temperature 


January 










May 












July 

September 

October 






























High 

Moderate 

Low 

Not enough data or not applicable 
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light, nutrient availability, and surface production. General year to year water mass characteristics outside 
the sanctuary are significant contributors to the overall relationship between variables and more study is 
needed to fully address local and regional variability of surface properties and their impacts on living marine 
resources. 

3.7 Interannual Patterns of Surface Properties as Indicators to Environmnental Perturbations 

Ocean climate anomalies for the PNW region reflect deviations from expected (normal) conditions 
that are sometimes strongly affected by interannual variability in climate. Through climatological anomalies, 
seasonal variability is, in part, removed, so that the anomalies can be assessed in relation to environmental 
perturbations with a higher degree of confidence. In this section, patterns of local to regional interannual 
variability of surface water properties will be presented through climatological anomalies using SST, SSHA, 
chlorophyll and turbidity data. General patterns involving anomaly patterns for these surface level properties 
will be visually interpreted for extremes (positive and negative) during the period of record of each data- 
set. Comparisons will be made with heightened periods of ENSO using the Multivariate ENSO Index from 
NOAA's Earth System Research Laboratory, Climate Analysis Branch (Wolter and Timlin, 1998). 

Interannual patterns of SST and SSHA 

SST Hovmoller plots reveal anomalous conditions inside the sanctuary during the well documented 
ENSO warm (Figure 38, orange) El Nino phases in winter 1997-1998 and 2003 and the cool (blue, Figure 
38) La Nina phases in 1988-89 and 1998-99. Cold events can be associated with positive ENSO index 
values (Figure 38, SST hovmoller), as was the case in July 2002. This period featured moderate to strong 
negative anomalies and an intrusion of cold subarctic waters (Grantham et al. 2004). SST anomaly maps 
for warm and cold ENSO and non-ENSO periods are depicted (Figure 39) with mainly positive (top row, 
orange patterns) anomalies evident during warm El Nino. During cool La Nina phases, negative anomalies 
are common (middle row, blue patterns, Figure 39). During non-ENSO periods, anomalous patterns tend to 
be slightly more spatially limited and shorter in duration. For example, extreme positive anomalies occurred 
in September 1990, August 1994, and negative anomalies were evident in December 1985 and June 1991. 
Typically, these patterns result from large-scale physical circulation fluctuations related to the CCS (Goericke 
etal. 2005). 
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Figure 38. Time-series plot (Hovmoller) along a North-South transect line revealing weekly gap-filled CoRTAD SST, and weekly 
SST anomaly values for the period 1985-2005. 
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Figure 39. CoRTAD monthly mean SST anomaly image maps during ENSO and non-ENSO periods. 



Sea surface height deviation was particularly evident during the ENSO warm (Figure 40, orange) 
phase in winter 1997-1998, and the strong anomalously cool (Figure 40, blue) spring of 1999. A general lag 
associated with SSHA temporal deviation from the end of El Nino through the start of La Nina is evident from 
the Hovmollerand the image maps, where heights remained positive until late 1998 (Figure 40). Spatial pat- 
terns of both negative and positive SSHA deviation revealed strong regional signals for most of the eastern 
North Pacific basin (Figure 41). Anomalously lower sea heights were evident from January 1999 until late 
summer 2002. Also evident is the moderate El Nino in 2003 (top row, Figure 41). Similar to SST anomaly, 
during non-ENSO periods, SSHA temporal deviation patterns tend to be slightly more spatially limited and 
shorter in duration (bottom row, Figure 41). 
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Figure 40. Time-series plot (Hovmoller) along a North-South transect line revealing (a) weekly AVISO SSHA; 
(b) climatological weekly SSHA; and (c) weekly SSHA temporal deviation for the period 1 992-2006. Multivari- 
ate ENSO index values added for reference. Source: NOAA/NWS. 
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Figure 41. AVISO monthly mean SSHA temporal deviaton image maps during ENSO and non-ENSO 
periods. 
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Interannual Patterns of chlorophyll and turbidity 

Monthly chlorophyll concentrations are strongly influenced by interannual variability in climate. The 
phase and strength of climate signals, such as the Pacific Decadal Oscillation (PDO), ENSO, and the North 
Pacific Gyre Oscillation (Thomas et al. 2003, Sackmann et al. 2004, Di Lorenzo et al. 2008), can lead to 
large-scale anomalies in chlorophyll concentration. This section presents patterns of local to regional interan- 
nual variability of climatological anomalies using SeaWiFS chlorophyll and turbidity data. General patterns 
involving stronger anomaly signals will be visually interpreted for select months in the 10-year climatology. 

Specific OCNMS anomalous chlorophyll "events" for the period were in the fall of 1999, September 
2000, and May and June 2007 (Figure 42). A strong negative chlorophyll anomaly was measured in the fall 
of 1997, associated with the strong El Nino event. Weaker negative anomalies were also evident in the fall 
of 1998, 2002, 2003, and 2005. Image maps generally representing warm and cool periods of ENSO are 
depicted in Figure 43, with strong positive anomalies generally associated with La Nina conditions, and nega- 
tive PDO. During warm ENSO phases, open ocean characteristics in the North Pacific basin reveal positive 
chlorophyll to the west of the study area, with negative anomalies inside the sanctuary. Strong non-ENSO 
events occurred in October 2001 and the summer of 2006, when anomalous chlorophyll was evident along 
the entire length of the study area (Figure 43, bottom row) that persisted through August 2006. Overall, 
chlorophyll anomalies appear to have only slight association with the large scale climate indices, such as the 
Multivariate ENSO Index (MEI) and PDO Index, particularly at localized scales. 
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Figure 42. Time series plot (Hovmoller) along a North-South transect line revealing filtered SeaWiFS chlorophyll and 
chlorophyll anomaly values for the period 1997-2007. Multivariate ENSO index values added for reference. Source: 
NOAA/NWS. 



An Oceanographic Characterization of the Olympic Coast National Marine Sanctuary 



ENSO- Warm Phase 

129°W 126°W 123°W 




Chlorophyll Anomaly (|wg/L) 







+0.5 



Figure 43. SeaWiFS monthly mean chlorophyll anomaly image maps during ENSO and non-ENSO periods. 



Monthly turbidity anomaly levels are extremely moderate for the region, but vary according to interac- 
tions with other physical drivers related to seasonal rainfall/runoff patterns. Higher values were observed 
in the 1998, 1999, 2001, 2006, and 2007 winters south of the Juan De Fuca Strait (Figure 44). Small-scale 
anomalies in Rrs670 were observed within the OCNMS in the winter-spring of 1999 and 2007 during La Nina 
conditions (Figure 45, middle row). Turbidity anomalies are mainly restricted to areas along the mid- to inner 
shelf and can span consecutive months. This was the case in early 1999, when elevated turbidity occurred 
over a 4-month interval inside the sanctuary (Figure 44). An abnormally low turbidity winter occurred in 
2005. 
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Figure 44. Time series plot (Hovmoller) along a North-South transect line revealing filtered SeaWiFS turbidity and 
turbidity anomaly values for the period 1997-2007. Multivariate ENSO index values added for reference. Source: 
NOAA/NWS. 



ENSO- Warm Phase 

129°W 126°W 123°W 




ENSO- Cool Phase 



z 


&S:*M ^5@B 


CO 






, 3^j 




■'■*?*=.■■. T"* ^M 




.j*- Mm 


Z 












^ $■ 


z 




c\J 


tP 



„-**f 


c ~"- %HBk9 


^< 





Non- ENSO- Positive Anomaly 





Turbidity (Rrs670) Anomaly (sr-1) 
-0.0005 +0.0005 +0.005 



Figure 45. SeaWiFS monthly mean turbidity anomaly image maps during ENSO and non-ENSO periods. 
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Summary of Surface Anomalies for the OCNMS 

Ocean climate anomalies developed from long-term data series were found to be associated with 
both large-scale climate cycles and localized perturbations in OCNMS environments. Selected data extrac- 
tions of climate anomalies inside the sanctuary are summarized (Table 4) to reveal general coincidence with 
ENSO indices. For example, during the strong El Nino period in fall 1997 through summer 1998, high SST 
and increased vertical SSHA deviation were detected inside the sanctuary for all months selected. Chlo- 
rophyll concentrations revealed strong negative anomalies at the beginning of the period, with alternating 
patterns at the end of the period. Turbidity anomalies revealed negative anomalies at the beginning of the 
period, with slight positive anomalies at the end. During the strong La Nina period in 1999, SST and SSHA 
were found to be negative, while chlorophyll and turbidity were strongly positive. During a 6-month period in 

2001 , negative anomalies were recorded for all variables analyzed. The ENSO index was positive again in 

2002, with unexpectedly colder water (negative SST anomaly) evident and decreased sea surface heights 
inside the sanctuary5. Lower turbidity and less chlorophyll production were evident inside the Sanctuary 
during the 2002 fall period. In 2004, neutral ENSO conditions were evident, with warmer than expected SST 
and SSHA, coupled with stronger than normal chlorophyll present in the sanctuary. Significant alternating 
patterns of ocean color variability reveal the potential volatility of ocean color signals inside the sanctuary. 
The extreme positive and negative anomalies of SSHA and SST are usually well separated and can last for 
years, generally revealing similar patterns to the ENSO index. However, heightened ENSO does not always 
imply extreme anomalous surface patterns, as seen in 2002 and 2004 (Table 4). 

Table 4. Comparison of ENSO index values and remotely sensed variable anomalies from extracted data averaged for the Olympic 
Coast National Marine Sanctuary 
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4. CONCLUSIONS AND RECOMMENDATIONS 

This work was intended to summarize surface layer conditions (both typical and atypical) as depicted 
from long-term satellite data series. The suite of remote sensing data, presented as conventional maps, 
Hovmoller diagrams, and statistical plots, afforded simplicity in summarizing conditions for the region. In 
particular, merged satellite oceanographic data revealed information on: 

• seasonal cycles and transition periods 

• linkages between seasonal warming and cooling, upwelling/downwelling, transport, and wind influences 

• seasonal nutrient sources and physical limiting factors controlling surface layer responses 

• interannual cycles and general response to climate oscillations (ENSO) 

All surface layer data collected and analyzed near the Olympic Coast National Marine Sanctuary 
exhibited change during a typical annual cycle. Using climatological summaries of average (expected) condi- 
tions, and through correlation maps and Hovmoller diagrams along three respective transect lines, periodic 
seasonal and subseasonal shifts in surface conditions affecting the OCNMS were interpreted. Five fairly dis- 
tinct seasonal cycles were evident in the ocean climate data summaries and were consistent with the known 
literature. 

Insights on underlying processes and limiting factors controlling surface response and nutrient sourc- 
es affecting the sanctuary were documented. The correlation patterns of surface water properties were found 
to be spatially and temporally tied to upwelling/downwelling cycles, outflow from the Juan de Fuca Strait, the 
Columbia River, coastal/estuarine influence, and increased light availability (Figure 37). 

There is general agreement that January is typically dominated by downwelling favorable winds, pre- 
cipitation influence (coastal/estuarine), and Columbia River discharge. During May, strengthened upwelling, 
increased regional heating and Juan de Fuca outflow influence surface layer characteristics most dramati- 
cally. In July, surface layer responses are indicative of enhanced upwelling, increased light and temperature, 
with available nutrients out of the Juan de Fuca Strait promoting chlorophyll production along the entire 
Olympic Coast. Late September/early October is typically dominated by increased chlorophyll production 
and westward propagation of SST fronts and eddies. Increased Ekman transport, driven by local alongshore 
winds, appears to be the governing forcing mechanism on this early fall chlorophyll peak enhancement 
(Strub, pers. comm.). Numerous surface reflections appear to be coincident with primary chlorophyll peak 
including decreases in localized SST, suppressed SSHA, and lower turbidity conditions that are evident in a 
2-week interval prior to and during chlorophyll peaks. Further investigation is needed to assess the periodic- 
ity and predictability of this phenomenon. 

The remote sensing data also provided insights on isolating upwelling responses and the nature of 
frontal features for the region. Since the surface patterns that were represented revealed periodicity and co- 
variability, further testing of relationships could be used to discern specific cues involving upwelling response, 
based on a combination of surface and subsurface characteristics. In addition, chlorophyll variability indica- 
tors were identified as a proxy for phytoplankton activity at specific time periods during the upwelling cycle. 
This information could also be used to discern specific timing cues for individual species assemblages. 

This report reviewed and summarized remotely sensed proxies needed for building appropriate deci- 
sion support tools, including integrated assessments. General uses of core variables estimated from satellite 
are shown (Figure 46) as a basis for environmental characterization, resource evaluation and problem detec- 
tion. Characterization and resource evaluation involve developing baselines of environmental parameters for 
trend detection, where data is extracted and subsequently integrated with biological and habitat information. 
Problem detection involves using remote sensing data to identify physical oceanographic and climatologi- 
cal conditions (and/or processes) directly and indirectly related to coastal water quality issues (e.g. low DO 
conditions and harmful algal blooms). It is suggested that a thorough analysis of expected surface conditions 
using a suite of merged information is critical in characterization and resource evaluation, problem detection, 
and trend determination. Assimilation of remotely sensed data for environmental modeling purposes will con- 
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tinue to evolve through improvements in satellite 
engineering and optics. 

An essential requirement for marine re- 
source evaluation and assessment involves the 
integration of spatially and temporally resolved 
oceanographic datasets. Through this work, a 
merged suite of oceanographic data was sum- 
marized to provide a foundation for building as- 
sessment capability and design, wherein merged 
data types are of more significance than singu- 
lar data types. Integration of this form works in 
many ways to improve remote sensing utilization 
through a more spatially and temporally compre- 
hensive assessment of surface patterns and vari- 
ability, and an increased understanding of interac- 
tions between processes, habitats, and resources 
necessary for sustained ecosystem function. De- 
velopment of indicators of biological response us- 
ing a full suite of remote data types should be a 
priority for future work involving marine resource 
evaluation and assessments. 
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Figure 46. Schematic describing general uses of remotely sensed 
parameters (proxies) outlined in this report for ecosystem resource 
evaluation, problem identification, and management decision sup- 
port. 



NOAA's National Center for Coastal Ocean 
Science (NCCOS), Office of National Marine 
Sanctuaries, and National Marine Fisheries Ser- 
vice have consistently worked to merge datasets (and disciplines) for improved ecosystem-based manage- 
ment. As an example, NCCOS's Biogeography Program relies on multifaceted data to study the relationship 
of species distribution patterns relative to geographical differences in the environment, with added specificity 
to marine-managed boundary concepts (NCCOS, 2005). For many protected areas, basic biological and 
physical data are lacking in spatial and/or temporal resolution. This remote sensing baseline is central for 
understanding the region in a broader context and how the biology and oceanographic conditions of a region 
change through time. 

In summary, the data and analysis presented here were intended to be the first of many steps to ad- 
dress present and future management questions posed by the OCNMS. The intent was to demonstrate how 
long-term data series could contribute to a better understanding of oceanographic variability and processes 
along the Pacific Northwest and elsewhere. This summary should be used as a vehicle for future discussions 
and work posed by joint partnerships between funding agencies, including NOAAand the research commu- 
nity vested in ocean science. 
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